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Abstract 

A full calculation of lepton-pair angular characteristics is carried out for e + e~ 
pairs created in pp, pn and pd collisions at intermediate energies. It is demonstrated 
that the proposed new observable, the dilepton decay anisotropy, quite sensitively 
changes for different sources and may be useful for their disentangling. The rele- 
vance of the dilepton decay anisotropy is shown in the context of a puzzling energy 
behavior for the ratio of the lepton yield from pd to pp reactions as observed at the 
BEVALAC. 
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As is known since a long time dileptons are quite attractive probes since they provide 
almost direct information on hot and dense nuclear matter during its evolution in heavy- 
ion collisions at BEVALAC/SIS and SPS energies |]J. The information carried out by 
leptons may tell us not only about the interaction dynamics of colliding nuclei, but also 
on properties of hadrons in the nuclear environment or on a possible phase transition of 
hadrons into a quark-gluon plasma. 

However, there are a lot of hadronic sources for dileptons because the electromagnetic 
field couples to all charges and magnetic moments. In particular, in hadron-hadron col- 
lisions, the e + e~ pairs are created due to the electromagnetic decay of time-like virtual 
photons which can result from the bremsstrahlung process or from the decay of bary- 
onic and mesonic resonances including the direct conversion of vector mesons into virtual 
photons in accordance with the vector dominance hypothesis. In the nuclear medium, 
the properties of these sources may be modified and it is thus very desirable to have ex- 
perimental observables which would allow to disentangle the various channels of dilepton 
production. 

For a decomposition of e + e~ sources contributing to the dilepton invariant mass spec- 
tra it seems to be natural to start from the study of elementary nucleon-nucleon collisions 
and then to move successively to nucleon-deuteron and more complicated systems. A step 
in this direction has been done in Ref. @ where the pd/pp dielectron ratio was measured 
in the 1 — 5 GeV energy range. In contrast to naive expectations, this pd/pp ratio displays 
a puzzling beam-energy dependence which is still a matter of debate. There are several 
interpretations of this effect: the interferences of different channels [0], the contribution 
of the inelasticity effects and of the subthreshold n -meson production ||. Even by 
considering the new and more precise data on the transverse momentum distributions of 
lepton pairs || the puzzle could not be resolved unambiguously. 

Recently, we have proposed to use lepton pair angular distributions for a distinction 
between different sources || [/]]. Indeed, the coupling of a virtual photon to hadrons 
induces a dynamical spin alignment of both the resonances and the virtual photons. One 
thus can expect that the angular distribution of a lepton will be anisotropic with respect 
to the direction of the dilepton (i.e. virtual photon) emission. It has been shown that due 
to the spin alignment of the virtual photon and the spins of colliding or decaying hadrons, 
this lepton decay anisotropy turns out to be quite sensitive to the specific production 
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channel §, % § . 

Following our previous work ||, |7j we focus in this letter on the study of the dilepton 
angular anisotropy as a way to disentangle various e + e~ sources in pN and pd reactions. 
Alongside with an extended computation of the Dalitz decay channel for pseudoscalar 
mesons, we present the first full results for nucleon-nucleon collisions at energies of a few 
GeV. By a simple extension of these calculations to the nucleon-deuteron case we will 
investigate if the anisotropy effect can help in solving the pd/pp dilepton ratio puzzle. 

As in |5], |7| we choose the polar 9 and azimuthal ip angles of the momentum /_ of 
a created electron with respect to the momentum q of a virtual photon to characterize 
the decay anisotropy, where l-,l+ are measured in the rest frame of this virtual photon, 
i.e. q = /_ + / + = 0. For comparing the shape of the angular distributions for different 
channels, the differential cross section for dilepton production in the channel i may be 
presented in the form: 



Si(M, 9) ee — p = Ml + * cos 2 9), (1) 
dm d cos 6 

where M is the invariant mass of a lepton pair (M 2 = g 2 — <f) and is defined by the 
normalization to unity of the angular distribution of ([!]) . The decay anisotropy coefficient 
for channel i, i.e. B^, then is given by 

_ ft(M,fl = 0°) _ 

~~ Si(M, 9 = 90°) [ ) 

Since the coefficient Bi is sensitive to the spin structure of the interacting hadrons, it is 
in general a function of M and the masses of the hadrons involved in the reaction. 

The total differential cross section for proton-nucleon collisions can be represented as 
a sum of differential cross sections for all channels: 

daPN = y ^ = A(l + BP N cos 2 9). (3) 

dMd cos 9 i= t nnel dMd cos 9 1 ; K> 

For the total anisotropy coefficient we then get 

do; Bi 



BPN= ^ <BfN>] <B f>= . *, » W 



i=channel ^ ^ 
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where the special weighting factors originate in the necessary angle-integrations. Thus, the 
anisotropy coefficient B pN is the sum of the "weighted" anisotropy coefficients (< Bi >) 
for each channel i obtained by means of the convolution of Bi with the corresponding 
invariant mass distribution. 

In a preceding study the anisotropy coefficients were calculated for the bremsstrah- 
lung and A-Dalitz decay channels within a microscopic One-Boson-Exchange (OBE) 
model J|. We will use the latter results in the present investigation without any modi- 
fication. Note that we take into account the interference terms between bremsstrahlung 
and A-channels for proton-nucleon and proton-proton interactions. As a result, the 
contributions of these two channels to (§) are reduced to a single coherent term. 

For the Dalitz decay of an 77-meson in Ref. @ a first simple estimate B n = 1 had been 
made appropriate for an ?7-meson with vanishing momentum in the center-of-mass system 
(cms) of the colliding hadrons. In this case the direction of a virtual photon q v is not 
influenced by the 77-meson momentum P v which simplifies the kinematical considerations. 
In extension to we now include the full dynamics, P v 7^ 0, where the direction of the 
dilepton momentum in the cms of the colliding nucleons, q, does not coincide with that 
in the eta rest frame q v ; as a consequence the anisotropy coefficient becomes a function 
of the invariant mass M and the collision energy T^. 

To take this functional dependence into account the //-channel for the lepton differen- 
tial cross section (|I|) is written in the following form: 

r -> da r i 

S v (M,e) = dP,-^ \T(M,9,P V )\ 2 , (5) 
J dP v 

where the averaged transition matrix element of the //-decay is given by 

\T(M, 9, P v )\ 2 ~ J dV -^j {e^miy e uapa q p qi a ) ■ L^, (6) 

and q, q\ are the four momenta of a virtual and real photon, respectively. The integral is 
taken over the available phase-space volume. The lepton tensor L^ u is defined as 

Lpn, = Tr /_7 At / + 7 !y , (7) 

where / = 7 CT / CT , /+ are the four momenta of the leptons. 

For the rj production cross section da v /dP v we use the phase-space oriented expression: 



with Ml = (p a + pb — P v ) 2 = s — 2y/sE v + nit. Here s — (p a + pb) 2 is the center- 
of-mass energy squared and p a ,pb,P v are the four- moment a of the colliding nucleons 
and the //-meson, respectively, while E v is the energy of the //-meson in the cms of 
the colliding nucleons. It is noteworthy that the particular case considered in Ref. P 
corresponds formally to the substitution da 71 /dP v ~ S(P V ). The modification of the phase- 
space formula (§) is included in the term T^(P V ), which is the production matrix element 
squared. Here we use a simple parametrization of the calculations presented in Ref. ||, 
where the //-production cross section in nucleon-nucleon collisions was calculated on the 
basis of an effective OBE model. 

Results of our numerical calculations for B v at initial energies from the rj production 
threshold to about 5 GeV are shown in Fig. [I]. We find that the anisotropy coefficient 
indeed depends on both the initial energy E and the dilepton mass M. B v approaches 
1 only close to threshold or for M — > 0. This behavior is quite natural because at 
the threshold energy we have \P V \ ~ 0; for low mass dileptons their velocity is close to 
the velocity of light and thus the difference between the two directions discussed above 
becomes negligible. If the invariant mass of a virtual photon is close to its kinematical 
limit (M — > m,), the anisotropy coefficient approaches a lower limit of about 0.1 -7- 0.25 
for bombarding energies from 1.26 to 5 GeV. 

In Fig. ||the dilepton invariant mass distributions da/dM from pp and pn interactions 
are shown at bombarding energies of 1.26 GeV and 2.1 GeV, which are the necessary 
ingredients for calculating the weighted anisotropy coefficients, Eq. (|9|). In Fig. || the 
"77" denotes the contribution of the ^-channel the "A" labels the contribution of the A 
resonance term, while u Br" denotes the bremsstrahlung channel. The dashed curves are 
the sum of bremsstrahlung and A-channel with interference. 

Using the pN cross sections for A-production and bremsstrahlung channels from 
Ref. P], we are now in the position to present the first estimate for the weighted co- 
efficients < E>i(M) > for pp and pn collisions at the energies of 1.26 and 2.1 GeV (Fig. |3]). 
The notation is the same as in Fig. [2|. The solid line denoted by "a//" represents the sum 
of bremsstrahlung, A-channel with interference and //-decay source. At the energy of 
1.26 GeV the anisotropy is determined by the A-decay since its contribution is dominat- 
ing at this energy (cf. Fig. |2|). The large contribution of the //-source at 2.1 GeV leads 
to the specific structure in < E>i(M) > for M < 0.3 GeV. 



5 



As mentioned above, a possible explanation for the experimental pd/pp dilepton ratio 
is related to the contribution of leptons from n-decay || . Whereas the energy threshold 
of ^-production in pN interactions is about 1.26 GeV, in the proton-deuteron case the 
effective threshold is lower due to the Fermi motion of the constituent nucleons in the 
deuteron. Thus at projectile energies below 1.26 GeV the contribution from the 77-Dalitz 
decay to the dilepton spectrum may be significant for the pd case, but negligible for pp. 
Such an enhancement of the 77-yield in pd compared to pp reactions has been also seen 
experimentally by the PINOT collaboration [ItJ. In line with this argument we expect 



larger values for the pd/pp dilepton ratio at low energy, and according to the calculations 
above one can expect that subthreshold 77-production may manifest itself also in the 
angular characteristics (anisotropy coefficients) of the created dileptons. 

In a first order approximation the decay anisotropy coefficient for pd collisions may 
be written as the sum of Bi for pp and pn interactions. However, this approximation is 
not sufficient in the energy region close to the ^-threshold. To take the contribution to 
the dilepton spectrum from near-threshold 77-production into account we calculate the 
coefficient B pd as 

B pd = £< Bg A > + < Bf >, (9) 

k=p,n 

da pd Bf 



dM 1 + y 

K K >= r 1 " ** 1 • (10) 

C^ 1 dM i + 1rw, dM i + Irp^ 

k=p,n J- \ ^-Ofir.A L ' g Dr i 

It should be stressed that in Eqs. (Q), (|IUD we take into account the interference between 
bremsstrahlung and A-channel for pp and pn separately. 

Fig. H shows the resulting weighted anisotropy coefficients < B,i{M) > for pd collisions 
at 1.26 and 2.1 GeV. The cross section da pd /dM is taken from Ref. f| and the Paris 
deuteron wave function |TTJ is used for describing the internal nucleon distribution in the 
deuteron. As expected, at threshold energy the total < B a u(M) > coefficient for the pd 
reaction differs substantially from that for pp or pn reactions due to the near threshold 
n-dec&y source. For higher energies (2.1 GeV) the shape of < B a u(M) > is similar in all 
these cases. This energy-dependent effect is more pronounced in Fig. |5| where the ratio 
of the anisotropy coefficients for pd to pp reactions, Rb = B pd /B pp , is represented at 
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the energies of 1.26 and 2.1 GeV. The comparison of the solid (including the r^-Dalitz 
decay) and dashed curves (without the ^-channel) illustrates the relative importance of 
the 77-decay contribution to the observable decay anisotropy of the dileptons. 

Thus, summarizing, the calculated anisotropy coefficients for pp, pn and pd reactions 
support our suggestion in Refs. |], [7] that the dilepton decay anisotropy may serve as an 
additional observable to discriminate the dilepton sources by experimental means. 

This research is partly performed in the framework of the Grants N°_ MP8300 and 
RFE300 from the International Science Foundation and Russian Government. The work 
of V.D.T was supported by Grant N° 3405 from INTAS (International Association for 
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Figure captions 



Figure 1: The anisotropy coefficient B v for the ^-channel at initial energies from 1.26 to 
4.9 GeV. 

Figure 2: The cross sections da/dM for pp and pn interactions at 1.26 and 2.1 GeV bom- 
barding energies. The "77" denotes the contribution of the ^-channel, the "A" labels the 
contribution of the A-resonance term, while Ll Br" denotes the bremsstrahlung channel. 
The dashed curves are the sum of bremsstrahlung and A-channel with interference. 

Figure 3: The weighted anisotropy coefficients < Bi(M) > for pp and pn collisions at 
1.26 and 2.1 GeV. The notation is the same as in Fig. [| The solid line denoted by 
"all" represents the sum of bremsstrahlung, A-channel with interference and 77-decay 
contributions. 

Figure 4: The weighted anisotropy coefficients < Bi(M) > for pd interactions at 1.26 and 
2.1 GeV. The notation is the same as in Fig. 0. 

Figure 5: The ratio of anisotropy coefficients for pd to pp reactions at energies of 1.26 
and 2.1 GeV. The solid lines correspond to calculations taking into account the 7/-decay 
contribution while the dashed lines do not include the ?7-Dalitz decay. 
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